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We derive high-resolution upper bounds for optimal product quantization of pathwise continuous 
Gaussian processes with respect to the supremum norm on [0,T] d . Moreover, we describe a 
product quantization design which attains this bound. This is achieved under very general 
assumptions on random series expansions of the process. It turns out that product quantization is 
asymptotically only slightly worse than optimal functional quantization. The results are applied 
to fractional Brownian sheets and the Ornstein-Uhlenbeck process. 
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1. Introduction 

In this paper, we investigate the functional quantization problem for pathwise continuous 
Gaussian processes A = (X t )tei,I = [0,T] d , where the path space E = C(I) is endowed 
with the supremum norm. For any real separable space {E, || • ||) and r € (0, oo), optimal 
quantization means the best approximation in L r E (¥) of a random vector X : (£l,A,F) — > 
E by random vectors X:(il,A, P) — > E taking finitely many values in E. If N E N, 
card(A(f2)) < N, then X is called A-quantization. This leads to the minimal level N- 
quantization error defined by 

e N . r {X, E) := inf{(E||A - X\\ r ) 1/r : AA-quantization of A}, (1.1) 

provided A € L r E (P). When E — R d , this problem is known as optimal vector quantization 
and has been extensively investigated since the early 1950s, with some applications to 
signal processing and transmission (see Gersho and Gray [11]) and to model-based clus- 
tering in statistics (see e.g., Tarpey [25]). Beyond these classical applications, optimal 
quantization has been used as a space discretization device to solve nonlinear problems, 
such as those arising in optimal stopping theory (American-style option pricing, reflected 
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BSDE, Bally and Pages [2]), nonlinear filtering (Pages and Pham [22]), forward-backward 
SDE (see Delarue and Menozzi [5]) and SPDE (see Gobet et al. [12]). The mathemati- 
cal foundations are treated in Graf and Luschgy [13]. Much attention has been paid to 
the infinite-dimensional case. This is the so-called functional quantization of stochastic 
processes: the aim is to quantize some processes viewed as random vectors taking values 
in their path spaces. Recently, a first application of functional quantization to statistical 
clustering of functional data has been investigated (see Tarpey and Kinateder [26] and 
Tarpey et al. [27]). The simplest application of functional quantization as a numerical 
method is to use it as an alternative to Monte Carlo simulation, using the quadrature 
formula 

E(F(X)) m E(F(X)) = F(a)F(X = a), where a = X(ti), 

for sufficiently regular functionals F : E — > R. If X is an L r -optimal iV-quantization and 
F is Lipschitz continuous, then the induced error is bounded by [F]u P eN,r(X, E), r > 1. 
Some numerical applications are being developed for the pricing of path-dependent op- 
tions (such as regular Asian options) in various models using E = L 2 ([0,T],dt) (Black 
and Scholes, Heston, see Pages and Printems [23], Wilbertz [29]). However, many im- 
portant functionals of processes, like those related to barrier options or to options on 
maximum, are only continuous with respect to the sup-norm on E = C([0,T]). 

Let us now describe what we will call the product quantization scheme. Let X be a 
centered valued Gaussian random vector. Let £i,£2>--- be i.i.d. Af(0, l)-distributed 
random variables and let (fj)j>i be a sequence in E such that £jfj converges a.s. 

in E and 

oo 

* = (1-2) 

3=1 

Let us call such a sequence admissible for X. For background on expansions for Gaussian 
random vectors, the reader is referred to Bogachev [4] and Ledoux and Talagrand [17]. 
One checks that (fj)j>i is admissible for X if and only if (fj)j>i is a normalized tight 
frame in the reproducing kernel Hilbcrt space (Cameron-Martin space) H = H x , that 
is, {fj,j >1}CH and Ej>i(fhhf H = \\h\\ 2 H for all h G H (sec Luschgy and Pages [21]). 
Then a sufficient (but not necessary) condition is that (fj)j>i is an orthonormal basis of 
H x . 

For m,N\,.. . ,N m € N with YYjLi Nj < N, let £j be an //-optimal N .-quantization 
for £j, that is, (E|^ — £j| r ) 1 / 7 * = ejv j . r (^ J R). An L r -product iV-quantization of X with 
respect to {fj)j>i is then defined by 



m 

X:=x im,..,N m) := j2£ jfj 



(1.3) 
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and the quantization error induced by X is 

(E\\X -X\\ r )^ r . 

Note that if a, = then the codebook a = X(£i) of X satisfies a = {J^JLi a ifi : a e 

]lJLi"i} and 

(Emm\\X - a\A 1/r <(E\\X-X\n 1/r . 

The minimal iVth product quantization error is then defined by 

e<$™ d) (X,E) := wi{(E\\X - X\\ r y/ r : (fj)j>i 6 E N admissible for X, 

XL r -product iV-quantization w.r.t. (fj)}- (1-4) 

Clearly, we have 

e N , r {X,E)<e%™ d) (X,E). (1.5) 

We address the issue of high- resolution product quantization in E = C(I) under the 
sup-norm, which concerns the performance of X = x( Nl ''"' Nm ' under a suitable choice 
of the marginal quantization levels Nj and the behaviour of e^ r ° d \x,C(I)) as N — > 
oo. For a broad class of Gaussian processes, we derive high-resolution upper estimates 
for e^ r r ° d) (X,C(/)). Furthermore, we describe a product quantization design X which 
attains this bound. Combining these estimates with precise high-resolution formulas for 
eN,r(X,C(I)) (see Dereich et al. [6], Dereich and Scheutzow [7], Graf et al. [14]), one may 
typically conclude that 

e^; d \x,C(I)) = O((\og\0gNye N AX,C(I))), 

for some suitable constant c > 0. This suggests that the asymptotic quality of product 
quantization, which is based on easy computations, is only slightly worse than optimal 
quantization. The optimality of this rate for product quantization rate remains open, 
although one may reasonably guess that it is optimal. 

The paper is organized as follows. In Section 2, we derive high-resolution upper esti- 
mates for e^ r ° d \x,C(I)) under very general assumptions on expansions. Section 3 con- 
tains a collection of examples, including fractional Brownian sheets, Ricmann-Liouvillc 
processes and the Ornstein-Uhlcnbeck process. 

It is convenient to use the symbols ~ and sa, where a n ~ b n means a n /b n — > 1 and a n ~ 
b n means a n = 0(b n ) and a n = Q(b n ). Throughout, all logarithms are natural logarithms 
and [x] denotes the integer part of the real number x. 
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2. High- resolution product quantization 

We investigate high-resolution product functional quantization of centered continuous 
Gaussian processes X = (X t )t^i on I = [0, T] d in the space E = C(J) equipped with the 
sup-norm ||a;|| = sup tg/ Let 

e^: A \x):=ef t f\x,C{I)). 

The subsequent setting comprises a broad class of processes. 
Let (fj)j>i £C(I) N satisfy the following assumptions: 

(Al) Wfj |j < Ci j-^ log(l + j)T for every j > 1 with > 1/2, 7 > and d < 00; 
(A2) fj is a-H61der-continuous and [fj] a < C2j h for every j > 1 with a <E (0, 1],6 £ R 
and C2 < 00, where 

(and |i| denotes the ^-norm of t £ R d ). 

In the sequel, finite constants depending only on the parameters T, 7, a, b, C±, C2, d 
and r are denoted by C and may differ from one formula to another one. Other depen- 
dencies are explicitly indicated. 

First, observe that by (Al), Y^jLi fj(t) 2 < X^/li ll/jll 2 < 00 f° r every t £ I, so wc can 
define a centered Gaussian process Y by Y t := X^=i £jfj(fi)- Using (Al) and (A2), we 
have, for p £ (0, 1], 

- fM = \Ms) - fM'M*) - fM 1 -" 
^dfjUs-tmWiW) 1 -" 

< C p j*'**'>-* 10g(l + j)T(l-P) \g - t \ aP 



and hence 



This yields 



3=1 



E|Y S - Y t | 2 = £ !/,(,) - / 3 (t)| 2 < ^[/^ |* - t| 2 ^ (2.1) 

3=1 \3'=1 / 

and using the Gaussian feature of Y, we obtain from the Kolmogorov criterion that Y 
has a continuous modification X. Consequently, (fj) is admissible for X and 

00 

3 = 1 
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For r £ [l,oo), the quantization error induced by the L r -product iV-quantization X 
X(N lt ...,N m ) ( sec ( L3 )) satisncs 



(E\\X-X\n 1/r = \\X-X\\ L r E{r) 

rn 

E& -&•)/:* 



< 



3=1 



- E 11& - &iU r (p)ii/j 
3=1 



E ^ 



j>m+l 



E ^ 



LJ(P) 



L-(P) 



so that 



(E\\X-X\n 1/r <Y,\\fi\\eN*AWO,l)) 
3=1 



E &/i 



j>m+l 



r\ 1/r 



For r £ (0, 1), we have 



(E\\X - X\\ r ) 1/r <E\\X- X\\ < E WfMSi - 01 +E 

3=1 

Let us now consider the truncation error. 



E fc/j 



j>m+l 



(2.3) 



(2.4) 



Theorem 1. Assume that (fj)j>i G C(/) N satisfies (Kl)-(K2). Then, for every n>2 
and r € (0, 00), 



E 



E^ 



3>" 



r\ 1/r 



< 



C(logn)T +1 /2 



(7 



#-1/2 



E 



E^ 



j">n 



r \ 1/r 



,#-1/2 ' 



Proof. By equivalence of Gaussian moments, 



E 



E^ 



j>n 



r\ 1/r 



< DE 



(2.5) 



for some constant D depending on r (cf. Ledoux and Talagrand [17], Corollary 3.2). The 
upper estimate for E|| J2j> n £j fj\\ ^ s based on corresponding estimates for finite blocks 
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of exponentially increasing length. For m > 1, set 

2 m 

Z = Z (m). = £ 

i=2 m - 1 +i 

For a given AT > 1, consider the grid Gat = {^j^- :i = 1, . . . , N} d . Then 
\\Z\\ < sup |Z t |+ sup \Z.-Z t \. 

teG N \s-t\KCN- 1 

It follows from the Gaussian maximal inequality that 



E sup \Z t \ <CJlog(l+N d ) sup yEZf. 
Using (Al), we have, for every t € I, 

2 m 2 m 

^z 2 t < e \\fj\\ 2 <c E r M iog(n-i) 27 <^2 m ( i - 2 *)m 2 T 

so that 



E sup \Z t \ <C0og(l+iV)2- m( ^ 1/2) m 7 . 
teGN 

Moreover, using (A2), we have, for |s — 1| < CW _1 , 



i=2"-i+l 

2 m 

< c| S -*r E 

J =2™- 1 +1 
2 m 

J =2™- 1 +1 

and hence 

2"* 

E sup \Z a -Zt\ <CN~ a E J 6 <CW~ a 2 m(1+6) . 
| 8 -t|<CJV-i 3=2—1+1 

Thus we have established the estimate 



E||Z( m )|| <C( v /log(l + A02- m (' 5 - 1 / 2 W + A^ a 2 m ( 1+fc )). (2.6) 
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As concerns the choice of N, set N := [2 um ] + 1, with u <G (0, oo) satisfying 1 + b — au < 
5—1?. Equation (2.6) then becomes 

Ej|Z (m )j| < C2- m ^- 1 '^m~< +1 l 2 . (2.7) 

We note that in the case b + 1} < —1/2, we may choose N = 1 and thereby obtain a 
power reduction from m 7+1 / 2 to m 7 . This can be improved. In fact, we have 

E\Zs-Z t \ 2 = \fj( s )-fj(t)\ 2 

2'" 

<C\s-t\ 2a j 2b <C\s-t\ 2a 2 m ( 1+2b ) 



so that 



d z (s,t) := (E\Z S - Z t \ 2 y/ 2 < C\s - t \ a 2 m{b+1 ' 2 \ 



If N{e,dz) denotes the covering numbers of I with respect to the intrinsic semi-metric 
dz, then, by chaining, 

rS , 

E sup \Z s -Z t \<E sup \Z s -Z t \<C y/logN{e,d z )de, 

\s-t\KCN- 1 dz{s,t)<8 JO 

where 6 := CN- a 2 m ( b+1 W (cf. Van dcr Waart and Wellner [28], page 101). Since 

'2"i(6+l/2) \ d / a 



N(e,d z )<C 



< £ < £0, 



and Jq 1 ^ylog(T/~x) dx < +00, we obtain, for sufficiently large N, 



^\ogN(e,d z )de <C2 m ^ b+1 ^ / ^/\og{l/x)dx <C2 m ^ b+1 ' 2 ^ 



Consequently, 



E||Z(" l )|| < C( Vlog(l + N)2- m ^- 1 ^m' 1 + 2 m ( fc+1 / 2 )) 

< C2 -m(,4-l/2) m -y ifb + -d<0. 



(2.8) 



We now complete the proof. For n > 2, choose m = m(n) > 1 such that 2 m 1 < n < 2 r 
Then 



3>n 



< E n^i 

j>m+l 



E^ 
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Since E|| E n <j<2"> O/j ll < E||^ (m) ll by the Anderson inequality (cf. Bogachev [4], Corol- 
lary 3.3.7), we deduce from equation (2.7) that 



E 



< 



j>m 



j 



,■7 + 1/2 



< 



Cm7 +i/2 ^ c(logn)^ +1 /2 



9.7(0-1/2) - 0™(^-l/2) 



,,0-1/2 



If 6 + < 0, then it follows from (2.8) that 



E 



j>n 



< 



C(logn) 7 

-0-1/2 ' 



Combining these estimates with (2.5) yields the assertion. 



□ 



Remarks. 

• The rate for the truncation error depends only on and 7, that is, on the decay of 
the size of functions fj (provided b + fl > 0). The occurrence of expansions with b + § < 
seems to be a rare event and otherwise b plays no role (see the subsequent example) . The 
case 7 = typically corresponds to one-parameter processes with / = [0, T] . 



The 



= (prod) 
'JV.1 



problem comprises the optimization of admissible sequences and, in 



view of (2.3) and (2.4), is thus related to the Z-numbers of X defined by 

r \ 1/r 



l n ,r(X) = l ntr (X,C(I)):=wi\ 



j>n 



(gj) admissible for X in C(I) 



(2.9) 



Rate-optimal solutions of the l„ ir -problem, in the sense of l rhr {X) ss (E|| Ej>n ^j9j\\ r ) 1 ^ r 
as n— > 00, have recently been investigated (see Kiihn and Lindc [16], Dzhaparidze and 
van Zanten [8, 9, 10], Ayache and Taqqu [1]). Admissible sequences of type (Al) and 
(A2) seem to be promising candidates. 



Example 1 (Weierstrass processes). Let 

fj (t) = j-* sm(j b +%, j > 1, > 1/2, beR,te [0, T]. 

Then \\fj\\ < j-® and [f^ = j b . Since / 3 -(0) = 0, we also have \\fj\\ < Tj b , so (Al) and 
(A2) arc satisfied, with 1? = max{i5, -6} and a = 1. The covariance function of X = 
: s;/; is given by 

3>1 
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Now, in the "Weierstrass case" b + $ > 0, we obtain, from Theorem 1, 

r\ 1 /r 



E 



j>n 



< 



,0-1/2 ' 



while in the "non- Weierstrass case," 6 + i9 < appears the better rate: 



E 



£&/> 



r\ 1/r 



< 



c 



-b-1/2 ' 



We pass to the minimal product quantization error e < ^ 1 ° d \x). 

Theorem 2. Assume that X admits an admissible set (fj)j>i in C(I) satisfying 
and ( A2 ). We then have, for every N > 3 and r £ (0, oo), 



(prod), y^gQgg log N)^ 7 



' N 



r(i)< 



(logiV)' 13 - 1 / 2 



(2.10) 



and 



(prod), C(l0g log 7V)^+7-V2 



'JV,r 



'P0< 



(logA^ 1 / 2 



«/ & + 1? < 0. 



Furthermore, the L r -product N -quantization X with respect to (fj), with tuning param- 
eters defined in (2.11) and (2.15) below, achieves these rates. 

Proof. Let r e [l,oo) and set Vj := j^ log(l + jo) 7 if j < jo '■= [e 7 ^] and Vj := 
j~' d log(l + j) 1 if j > jo- The sequence {vj)j is then decreasing. Since 

lim fcefe. r (7V(0, 1),R) exists in (0,oo) 

k — >oo 

(cf. Graf and Luschgy [13]), we deduce from (2.3), (Al) and Theorem 1 the estimate 



\j=l 



_! log(l + m) 7+1 / 2 



for every m,Ni,..., N m £ N with YYjLi Nj < -/V. (The case b + d < is treated analo- 
gously.) Consequently, 



e%?\x)<Ci*i 

li=i 



! log(l + m)T +1 / 2 



,0-1/2 



m,4..,JV m eN, 
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For a given TV e N, we may first optimize the integer bit allocation given by the Nj's for 
fixed m and then optimize m. To this end, note that the continuous allocation problem 
reads 



inf £ Vj yf : Vj > 0, f[ Vj < N \ = £ u.zj 1 = AT V™ m f[ Vj 



l/m 



o'=i 



where 



-l/m 



and Zi > • • • > z m . One can produce an (approximate) integer solution by setting 

/to \ 



(2.12) 



provided z m > 1 . Then 



l/m 



3=1 



Since the constraint on m reads m £ I(N) with 



-l/m 



/(TV) := \ m G N : iV 1 /" 1 ^ JJ ^ > 1 



we arrive at 

(prod) 



°>(X)<C inf 



mel(N) 



N- 1 / m \og(l + m)'f log(l + m)^ +1 / 2 



10-1 



,0-1/2 



(2.13) 



(2.14) 



for every N G N. We check that 7(iV) is finite, /(AT) = {1, . . . , m*(iV)}> m*(iV) increases 
to infinity and 

loe N 

m*(N) %- asN-*oo. (2.15) 



Finally, let 



such that 



2 log AT 

m = m(N)El(N), with m(N) < for N > 3 (2.16) 

log log N 



, \ 21ogA^ 

m(N) ~ ; — : — asA^^oo. 

v ; log log AT 
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This is possible in view of (2.14). Using (2.4), the case r € (0, 1) follows from r = 1 since 
the L r -optimal N 3 .-quantizations £j satisfy E|£j — £j| < CN^ 1 , j > 1; sec Graf et al. 
[15]. □ 

We may reasonably conjecture that for many specific processes, the above rate is 
the true one. This would imply that product quantization achieves the optimal rate for 
quantization, namely the rate of convergence to zero of eN,r(X) := eN >r (X,C(I)), only 
up to a loglog N term in formula (2.16). This is in contrast to the Hilbert space setting, 
where the optimal rate is attained by product quantization (cf. Luschgy and Pages [20]). 
To be precise, we summarize the results on eN,r{X) in the present setting. 

Proposition 1. (a) Assume that X admits an admissible sequence in C(I) satisfying 
(Al) and (A2). Then 

(loglog AT+V^ 
e N , r (X) = 0{ (logjvr , 1/2 ) (2-17) 

and 

e ^ (X) = Q ( w^ 2 )' <f* + *<°- (2.18) 

(b) Assume that X admits an admissible sequence satisfying (Al). Let /i be a finite 
Borel measure on I and let V :C(I) —> L 2 (I,ii) denote the natural embedding. Then 



e N AV{X), L 2 ( M )) = O 

and 

eff\v{X\L\»))=0 



(loglogiV)T 

(logA^-Va 
(log log N)i 



Proof, (a) The proof is not constructive. We use Proposition 4.1 in Li and Linde [18], 
which relates /-numbers (see (2.9)) and small ball probabilities (but this relation is not 
always sharp). By combining this relation and Theorem 1, we obtain 

/ / (<y+l/2)/(0-l/2)- 

- log(P(||X|| <e)) = 0[ s- 1 ^- 1 ^ log - 



/ / /iwT/C'-VaX 

•log(P(||X|| <e)) =0 e- 1/( ^ 1/2) logf - , iffe + i9<0 



as e — > 0. We may then apply a known, precise relationship between these probabilities 
and eN,r(X) (cf. Dereich et al. [6], Graf et al. [14]) and this leads to the desired estimate. 
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(b) Let (fj)j>i be an admissible sequence in C(I) for X satisfying (Al) and consider 
an L 2 -product TV-quantization of V(X) based on {Vfj)j>\, 



v(x) = E^(&). 

where £j are L 2 -optimal Voronoi TVj-quantizers; see Luschgy and Pages [19]. Then, using 
the independence of £j — £j, j > 1, and the stationarity property £j = E(£j|£j) of the 
quantization , we have 



E 



L*( M ) 



=E E l^-^l 2 ii^iii»(M)+ E iraiW) 

i=l j>m+l 

<c(EA7 2 i- 2 "io g (i+i)^+ ^ r M iog(i+i) 2 ^ 

\J=1 j>m+l / 

We then argue along the lines of Luschgy and Pages [19] to conclude that 



□ 



Sometimes, (2.17) provides the true rate for eN,r{X) (as for the two-parameter Brow- 
nian sheet), sometimes it yields the best known upper bound (as for the ci-parameter 
Brownian sheet with d > 3) and sometimes (2.18) provides the true rate (as for Brow- 
nian motion). The latter typically occurs when the rate of eM, r {X) and the "Hilbert 
rate" of eN^ r (V(X), L 2 (dt)) coincide (see Section 3). It remains an open question to find 
conditions for this to happen. 



3. Examples 



3.1. Fractional Brownian motions and fractional Brownian sheets 



We consider the Dzaparidze-van Zanten expansion of the fractional Brownian motion 
X = (Xt)te\o,T] with Hurst index p S (0, 1) and covariance function 



EX s X t = \{s 2p + t 2p -\s- t\ 2p ). 
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These authors discovered, in Dzhaparidze and van Zanten [9], that the sequence 
tiM TP Cp V2 . ( X jt\ 

U (*) = 77 — ; \, p+i sm r ' J - ls 



^^i^-ntjj. JS1 ' (3 - 1> 

in C([0,T]) is admissible for X, where J v denotes the Bessel function of the first kind 
of order < x\ < x^ < ■ ■ ■ are the positive zeros of J— p , < yi < y% < • • • the positive 
zeros of Ji_ p and c 2 p = T(l + 2p) wa.{np)/n. 
Using the asymptotic properties 

Xj ~yj~Tcj, Ji-pixj)^ J-p{yj) r 1/2 asj^oo 

7T 

(cf. Dzhaparidze and van Zanten [9]), one observes that a suitable arrangement of the 
functions (3.1) (like f%j = fj, faj-i = f?) satisfies (Al) and (A2) with parameters ■& = 
p + 1/2, 7 = 0, a = 1 and 6=1/2 — /?. Consequently, 

while (see Dereich and Scheutzow [7], Graf et al. [14]) 



e N ,r( FBM ) « (logN)-P. (3.3) 

The tensor products of functions (3.1) are admissible for the fractional Brownian sheet 
X over [0,T] d with covariance function 

d 
i=l 

Pi £ (0,1), and satisfy conditions (Al) and (A2) with d = p + 1/2, p = rrrmi<j<d/?j, 
7 = d{m — 1), where m = card{i G {1, . . . , d} : p% = p}, a = 1 and b = maxi<j<d(l/2 — 
Pi))+. This is a consequence of the following lemma which ensures stability of conditions 
(Al) and (A2) under tensor products. 



Lemma 1. For i e {1, . . . , d}, let (fj)j>i eC([0,T]) N satisfy (Al) and (A2) with pa- 
rameters ??i,7i,Oj,&i such that ji = 0. Then a decreasing arrangement of ((^) i=1 /j 4 )jeN ii 
satisfies (Al) and (A2) with parameters § = mini<j<d 7 = 7?(m — 1), where m = 
card{i S {1, . . . , d} : i?i = §}, a = mmi<j<d<ii and b = (maxi<i<d&i)+. 
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Proof. For j = (ji, ■ ■ ■ , 3d) € N d , set /, = ®ti f)< *° that fjjt) = nj=i * € 

[0,T] d . We have 

<2 d 

ll/,||<nil4H< c n^ and l/i(s)-/ j (i)|<C'maxjtk-r- 

— J — — l<2<d 
i=l i=l 

Let Wj := rii=i h f>i ■ Choose a bijective map -0 : N — > N d such that life := iw/b) is decreas- 
ing in k > 1. Set / fc := /,/,(£). Then 

u fe w Ck"^ (log fc)^'™- 1 ) as fc -> oo 

(cf. Papageorgiou and Wasilkowski [24], Theorem 2.1). Consequently, 

\\h\\<Ck-\\ogk)^ m -^ 

and, for j = ip(k), 

d d 



3i <J\ji< IK ^ ^(logfc)"^- 1 ' < Ck, 



»=i t=i 



hence 



i/ fc (s)-/ fc (t)i<cfc 6 | S -*r. D 



Therefore, by Theorem 2 and Proposition 1, 



3jv>r (FBS)-C^ ) (3.4) 



and 



e N , r (FBS) = O { j . (3.5) 

The Hilbert space setting E — L 2 ([0,T] d ,dt) provides the lower estimate 

e^FBS) = n ^ j (3.6) 

(see Luschgy and Pages [19, 20]). The true rate of eN.r(FBS) is known only for the case 
m = 1, where the true rate is the "Hilbert rate" (2.6) (see Dereich et al. [6]), and for the 
case m = 2, where (3.5) is the true rate (see Belinsky and Linde [3], Graf et al. [14]). A 
reasonable conjecture is that (3.5) is also the true rate for m > 3. 
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3.2. Riemann Liouville and other moving average processes 

For ijj <E L 2 ([0,T],dt) and a standard Brownian motion W, let 

X t = ( i>(t-s)dW s , t€[0,T], 



and assume that X has a pathwise continuous modification. Since 

ps/\t 

EX s X t = ip(s-u)ip(t-u)du, 



/ j W = V / |/*(*-»)c«(=<i^)d» 

i/w- Ctt-'f-' )^ i>l, (3.7, 

is an admissible sequence for X. Observe that (3.7) provides well-defined continuous 
functions, even for ift £ L 1 ([0,T],dt). 

Lemma 2. Let ip € L 1 ([0,1], dt). 

(a) Ifip(t) = Jq \ip(s)\ ds is (3-Holder continuous with (3 S (0, 1], then the sequence (fj) 
from (3.7) satisfies (A2) with a = f3 and 6=1. In particular, if ip E L 2 ([0,T],dt) , then 
( A2 ) is satisfied with a = 1/2 and b = 1 . 

(b) // "0 ftas finite variation over [0, T], then ( Al j is satisfied with $ = 1 and 7 = 0. 

Proof. Let Aj = (7t(j - l/2)/T)" 2 . (a) For s < t, we have 

= \/f { rVW(cos(( S - U )/VA7)-cos((t-u)/ v / Ai))d U 



f/^lt) COs((i — u)/ y'Xj) du 

J s 

so that 

1/3 (*) - / 3 -(*)i < yf (^iiv-iundt) + [ \m\ d«) ■ 

(b) We have 

/j (t) = - V /2A~/T jf d(sin( (t - s) / y/Xj) ) 

= y / 2VT^(0) sm(t/^X]) + J sin((f - «)/ ^) cty(s)) 
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so that 
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WfjW < y^/TdvKo)! + Var(V>, [0,T])). 
This lemma yields a universal upper bound, 



□ 



^(prod) 
-N,r 



(X) = 



log log N 
(log 1/2 



for functions ip having finite variation. 

In the sequel, we do not concern ourselves with improvements of the parameter b in 
(A2) since the condition b + $ < cannot be achieved in this setting. 

Lemma 3. Let ip e ^([O^df). 

(a) If ip is positive and decreasing on (0,T] and ip(t) = f*ip(s)ds is (3-Holder contin- 
uous with (3 £ (0, 1], then the sequence (fj) from (3.7) satisfies \\fj\\ < Cj~@. If [3 > 1/2, 
then (Al) is satisfied with d = (3 and 7 = 0. 

(b) // ip(0) = 0, ip is (3-Holder continuous with (3 € (0, 1] and ip is differentiable on 
(0, T] such that ip' is positive and decreasing on (0, T], then (Al) is satisfied with •& = 1 + (3 
and 7 = 0. 



Proof. Let Xj = (n(j - l/2)/T)" 2 . (a) For t < y/X], we have 

\f 3 (t)\< y/2/T<p(J\~). 
Using the second integral mean value formula, we obtain, for t 6 [y/Xj, T] and some 



\fj(t)\ < VVT 



ip(s)cos((t- s)/y/Xj)ds 



ip(s) cos((t — s)/y Xj) ds 



( fv^ /— 

2/Tl J ip{s)cos((t-s)/ v ^)ds 



cos((t — s)/y Xj) ds 



< V / 2/r(^( v /A,) + 2 v /A^-( v /A J )) 
<3^/2/T V (^X~). 
Consequently, 

||£||<3 v ^>(VA i )<C'A? 
(6) The function ^ is absolutely continuous on [0,T], so an integration by parts yields 



,/3/2 



f j (t) = y/2\ j /T / W)Bin((f- fl )/v/Aj) 



ds. 
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Arguing as in (a) (with -0 replaced by "0'), we deduce that 



(l+/3)/2 

3 ' U 



11/,-H <3 v /2a,/t^(Va7)<ca 

Now, let 0(t) =^" 1 / 2 with p £ (0,oo). Then 

X t =X£= [ (t-sy-^dWs, te[0,T] (3.8) 
Jo 

so that A p is a Riemann-Liouville process of order p. Using the (pA i)-H61dcr continuity 
of the application 1 1— > X p from [0,T] into L 2 (P) and the Kolmogorov criterion, we can 
check that X p has a pathwise continuous modification. 

Lemma 4. Let ?/>(£) = i p ~ 1 / 2 , p € (0,oo). TTien i/ie sequence (fj) from (3.7) satisfies 
(A2) with a = min{l,p + 1/2}, b=l and (Al) for p G (0,3/2] with d = p+ 1/2 and 
7 = 0. 

Proof. This is an immediate consequence of Lemmas 2 and 3. □ 
We deduce, for Riemann-Liouville processes of order p £ (0,3/2], that 

(prod) , „ M _ ^ / (^g A^)P+ V 2 \ 

ejVir (RL)-0^ (bgAOp j, (3.9) 

while for every p € (0,oo) (see [18], Graf et al. [14]), 

e N ^RL)^(logN)-P. (3.10) 

To go beyond p = 3/2, we must slightly change the way we quantize. Let rp(t) = 
tP- 1 / 2 , with p > 3/2, and choose k <E N such that k + 1/2 < p < k + 3/2. Set Xj = (7t(j - 
l/2)/T)~ 2 . For /c G {2n — l,2n} n £ N, integration by parts yields the expansion 



/iW = E(- i r" lA rV27T^ (2m " 1) W + (-l)"A; i v / 27r / ^ (2n) ( S )cos((i- S )/yA~)d S 

m=l Jo 

=: 9j {t)+hj{t), ie[0,T]. 

Since V (2n) (i) = C^" 1 if fc = 2n - 1 and V<» (t) = Ct 13 \ik = 2n with /3 = p-k-l/2e 
(0, 1] , we deduce from Lemma 2 and Lemma 3 that the sequence (hj) in C([0, T]) satisfies 
(Al) with •& = p+ 1/2, 7 = and (A2) with a = p - k - 1/2, & = -k if A; = 2n - 1 and 
a = 1, b = —k + 1 if k = 2n. Clearly, the sequence (gj) also satisfies the conditions (Al) 
and (A2) (with d = 2, 7 = 0, b = -2 and a = p - k - 1/2 if k = 2n - 1 and a = 1 if k = 2n). 
Consequently, there exist centered continuous Gaussian processes U = (Ut)te[o,T] and Z 
such that U = J2fLi£j9j a - s -> ^ = Y^jLi£.jhj a - s -> 

X = X P = U + Z (3.11) 
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and U € spanj-i^ 2 " 1-1 ) : m = 1, . . . , n} a.s. Observe that 

n 
m—1 

where Vm = J^f =1 A™^ is ^(O.E^i A 2m )-distributed. 

Now use, for example, [i\T 1 ' 2n ]-quantizations of rj m and a [viV]-product quantization 
of Z for the quantization of X (which is clearly not optimal in practise, but remains rate 
optimal). Let rj m be an i r -optimal [iV 1 ' 2 ™] -quantization for ?y m , 

n 

:= (-l) m " 1 V2/^V' (2TO " 1) ^m, 

m— 1 

and let Z^ be the //-product [V~N] -quantization of Z from Theorem 2. A (modified) 
//-product TV-quantization of X with respect to (fj) is then defined by 

X := + Z^. (3.12) 
Using Theorem 2, we can show for the quantization error, that 

\\U + Z-X\\ L ,<C(\\U-U^\\ Lr +\\Z-Z^\\ Lr ) 



<qX VVTU {2m - 1} \\\\v m - fuhr + \\z-z- 

C C(\og\ogVN) p+1/2 



Nu 



< 



< 



TV 1 / 2 ™ " (log VN)p 

C{log\ogN)P +1 / 2 
(log AT) p 



so that, with the above modification, (3.9) remains true for p > 3/2. 

Now, consider the stationary Ornstein-Uhlenbeck process as the solution of the 
Langevin equation 

dX t = -pXtdt + o-dWt, t£[0,T], 

with Xq independent of W and Af(0, ^ )-distributed, a > 0, f3 > 0. It admits the explicit 
representation 

X t = e- 0t Xo + ae-P* f e 0s AW s (3.13) 

and 



o 



EX s X t = ^e-^. 
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By Lemma 2, the admissible sequence 

satisfies conditions (Al) and (A2) with «?= 1, 7 = 0, a = 1 and 6=1. Consequently, 

( P rod), nm n ( log log AT \ 

(0 ^ )=0 i(bgjvF 5 J' (3 - 14) 

while (see Graf et al. [14]) 

e N , r {OU)™{\ gN)- l > 2 . (3.15) 
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